Spin-valley blockade in carbon nanotube double quantum dots 
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We present a theoretical study ol the Pauli or spin-valley blockade for double quantum dots in 
semiconducting carbon nanotubes. In our model, we take into account the following characteristic 
features of carbon nanotubes: (i) fourfold (spin and valley) degeneracy of the quantum-dot levels 
(ii) the intrinsic spin-orbit interaction which is enhanced by the tube curvature, and (iii) valley 
mixing due to short-range disorder, i.e., substitutional atoms, adatoms, etc. We find that the spin- 
valley blockade can be lifted in the presence of short-range disorder, which induces two independent 
random (in magnitude and direction) valley-Zeeman fields in the two dots, and hence acts similarly 
to hyperfine interaction in conventional semiconductor quantum dots. In the case of strong spin-orbit 
interaction, we identify a parameter regime where the current as the function of an applied axial 
magnetic field shows a zero-field dip with a width controlled by the interdot tunneling amplitude, 
in agreement with recent experiments. 

PACS numbers: 73.63.Kv, 73.63.Fg, 73.23.Hk, 71.70.Ej 
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I. INTRODUCTION 

Recent developments of experimental techniques al- 
low for preparation, manipulation and readout of few- 
electron spin states in quantum dots (QDs), 1 indicat- 
ing the strong potential of these systems for future ap- 
plication in quantum information processing!^ A ma- 
jor factor limiting the performance of quantum-dot spin 
qubits in widely used III-V semiconductors (e.g., GaAs) 
is spin decoherence due to hyperfine interaction with 
nuclear spins. A strategy to suppress spin decoher- 
ence is to use QDs dominantly consisting of nuclear- 
spin-free isotopes of group IV materials. Carbon struc- 
tures, such as carbon nanotubes (CNTs) or graphenc, 
are prime candidates for that purpose as the natu- 
ral abundance of spin-carrying 13 C nuclei is very small 
(1%). This observation has motivated intensive theoret- 
ical invest igatiorPEBl anc [ the experimental realization of 
QDs in carbon nanostructures.^"221 Further perspectives 
of carbon-based quantum information processing have 
been opened by proposals suggesting to utilize the val- 
ley d egree of freedom of the delocalized electrons as a 
qubitj22EU and to exploit the interplay of spin-orbit in- 
teraction, valley mixing, and the bending of CNTs for 
implementing qubit operations.^ 

The Pauli blockade or spin blockade effectP^ in con- 
ventional semiconductor double QDs (DQDs) has pro- 
vided a distinct probe of spin physics in these devices 
and has been utilized in the past decade for various pur- 
poses in the context of spin qubits. A basic application is 
spin-state initialization and readout in experiments real- 
izing resonant manipulation of single spins ! 34 " 36 ! Pulsed- 
gate techniques combined with the spin blockade setup 
have been usecpZH22l m qubit manipulation experiments 
where the information was encoded in the two-electron 
spin states S and To or S and T + . Similar experiments 
have been utilized to prepare the state of the nuclear-spin 
ensemble of the crystal lattice, with the aim of prolonging 
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FIG. 1: (Color online) Schematic of the spin-valley blockade 
setup with a carbon nanotube double quantum dot and an 
external magnetic field B aligned with the tube axis. In this 
regime electrons are transported from source (S) to drain (D) 
while the DQD occupancy changes between single and double. 
Spots represent electrons; the figure shows the (0,1) charge 
configuration of the double dot. Lead-dot tunneling rates Fl, 
Tr and interdot tunneling amplitude t are indicated. 



the decoherence time of the qubit HQJSS1 Furthermore, spin 
blockade has been proven an efficient tool to gain infor- 
mation about the mechanisms of spin relaxation and de- 
coherence, and the corresponding energy scales. In par- 
ticular, it has be en applied to measure the energ y scales 
of hyperfin e! 43 * 44 ^ and spin-orbit interactions! 45 * 46 * The im- 
plementation of this range of functionalities in carbon- 
based quantum dots, potentially showing improved qubit 
performance, is an intense ongoing effort! 20 " 22 * 29 * 

In this work , we consider Pauli blockade in a trans- 
port setupP^S where electrons are transmitted from the 
source to the drain in a serially coupled DQD via the 
(0,1) -> (1,1) -► (0,2) -► (0,1) cycle (Fig. [j). Here 
(nr,,nji) denotes the charge state with nr, (jir) elec- 
trons in the left (right) QD. In conventional semicon- 
ductor DQDs, if the (1,1) and (0,2) states are aligned 
in energy, then states sharing the same spin state be- 
come hybridized due to interdot tunneling. The only en- 
ergetically available (0,2) state has a singlet spin state, 
therefore it hybridizes with the (1,1) singlet only, leaving 



the three (1,1) triplet states without a (0,2) component. 
This implies that whenever a (1,1) triplet state is occu- 
pied in the transport process the current is blocked since 
the (1,1) triplet state cannot decay to a (0,1) state: the 
occupation of (1,0) states is energetically forbidden and 
the source connected to dot L cannot absorb the electron. 
The blockade can be lifted by various mechanisms influ- 
encing spin dynamics, e.g., hyperfme interaction, spin re- 
laxation, etc., resulting in a nonzero steady-state current, 
termed as the "leakage current" P through the DQD. 

This explanation is altered in the case of electrostati- 
cally defined CNT DQDs (Fig. [IJ where the val ley de gen- 
eracy of the electronic spectrum is maintaine d 10 ! 21 ! , re- 
sulting in fourfold-degenerate (spin and valley) quantum- 
dot energy levels (see Fig. [2]). In this case the six (0,2) 
states have combined spin-valley wave functions which 
are antisymmetric under particle exchange. We refer to 
such states as supersinglets. The 16-dimensional (1,1) 
subspace can be separated to a six-dimensional supers- 
inglet subspace and a ten-dimensional subspace of com- 
bined spin-valley wave functions being symmetric under 
particle exchange, i.e., supertriplets. Hybridization oc- 
curs between the (0,2) states and the (1,1) supersinglets, 
and the (1,1) supertriplets do not acquire any (0,2) com- 
ponents. As discussed above for the case of conventional 
semiconductor DQDs, this leads to a blockade of the 
transport. The blockade can be lifted by various mecha- 
nisms affecting the spin and valley dynamics. To distin- 
guish the cases of conventional and CNT DQDs, we refer 
to them as spin blockade and spin-valley blockade in the 
following, respectively. 

Here we focus on recent experiment ^2U22J observing the 
spin-valley blockade in clean CNT DQDs with natural 
(1%) and enriched (99%) 13 C abundance. Charge sensing 
data^ indicates that in these samples the valley degener- 
acy was maintained in contrast to other observations. 20 29 
In the case of the isotope-enriched samples, a zero-field 
peak has been observed^ in the magnetic-field depen- 
dence of the leakage current at small interdot tunnel- 
ing. Following a model developed for GaAs DQDspQ 
this feature has been attributed to hyperfme interaction, 
although the corresponding energy scale inferred from 
the measurement is two o rder s of magnitude larger than 
the theoretical estimates P^ This discrepancy has not 
yet been explained. At large interdot tunneling, th e ob- 
served magnetotransport data show a zero-field di p 1 21 * 22 ! 
with a width controlled by the transparency of the in- 
terdot tunneling barrier, irrespective of the dominant 
isotope species. In InA s DQ Ds, a similar feature has 
been measured recentlypSES and good agreement has 
been found with a phenomenological model incorporating 
spin-orbit-enabled spin-flip interdot tunneling and spin 
relaxation!^ The same mechanism might be responsible 
for the observed magnetotransport in CNT DQDs as well. 

In this work, we provide an alternative explanation of 
the zero-field dip found in the magnetotransport curve 
of CNT DQDs in the spin-valley blockade regime. Using 
a microscopic model, we argue that the disorder-induced 



valley dynamics is different in the two dots, resulting in 
the lifting of the spin-valley blockade and allowing for 
a finite current through the DQD. In this mechanism, 
disorder plays a role analogous to hyperfme interaction 
in conventional semiconductor DQDs. 

To show this, we set up a model Hamiltonian for the 
DQD accounting for the following unconventional prop- 
erties of CNT QDs: (i) fourfold (spin and valley) degen- 
eracy of the QD energy levels, (ii) the intrinsic spin-orbit 
interaction which is enhanced by the tube curvature and 
induces an energy splitting between Kramers pairs, and 
(iii) valley mixing due to short-range disorder, i.e., sub- 
stitutional atoms, adatoms, etc. We provide a micro- 
scopic analysis of property (iii), resulting in an effective 
Hamiltonian for a single fourfold degenerate QD level. 
We find that disorder appears in this Hamiltonian as a 
random (in magnitude and direction) effective magnetic 
field acting on the valley degree of freedom. We express 
this valley- Zeeman field as a function of the disorder con- 
figuration and the envelope function of the electron occu- 
pying the QD. Our transport calculations are based on a 
Born-Markov master equation. The main finding of this 
work is that the disorder-induced valley- Zeeman fields 
provide a mechanism that lifts the spin-valley blockade. 
Depending on the relative significance of spin-orbit inter- 
action and disorder, we identify different patterns in the 
magnetic-field dependence of the steady-state current. In 
the case of strong spin-orbit interaction, we find a zero- 
field dip in the magnetotransport curve, in agreement 
with recent experiments, however our model does not in- 
clude spin-flip tunneling or spin-relaxation processes. In 
the case of strong disorder, we find that the magneto- 
transport curve can show both a zero-field dip and peak, 
depending on the disorder configuration. 

The rest of the paper is organized as follows. In Sec. 
p"H we provide a microscopic analysis of short-range dis- 



order in a CNT QD. In Sec. Ill the model Hamiltonian 



of the CNT DQD and the master-equation approach is 
described. In Sees. [TV] and [Vj we study the magneto- 
transport in the spin- valley blockade regime in the cases 
of strong spin-orbit interaction and strong disorder, re- 
spectively. Our conclusions are presented in Sec. |VI| 



II. 



SHORT-RANGE DISORDER IN THE 
QUANTUM DOT HAMILTONIAN 



In this section, we consider a single electrostatically 
defined QD in a semiconducting CNT, in the presence of 
a homogeneous magnetic field. Our aim is to derive a 4x4 
effective Hamiltonian describing the effect of the short- 
range disorder present in the CNT on a fourfold (spin and 
valley) degenerate state of the QD. We show that in this 
effective Hamiltonian the short-range disorder appears as 
an effective magnetic (Zeeman) field acting on the valley 
degree of freedom, and having a random magnitude and 
direction. 

First we consider a CNT QD model without spin-orbit 
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FIG. 2: (Color online) (a) Schematic of a quantum dot in a 
carbon nanotube, with an external magnetic field B aligned 
with the tube axis, (b) Magnetic-field dependence of the spin- 
orbit-split single-electron ground state sublevels of a nanotube 
quantum dot, obtained from diagonalizing Hq + H c s, dis [see 
Eqs. H and pj]. Spin and valley quantum numbers of the 
energy levels are indicated on the right. 



interaction and short-range disorder. We choose the z 
axis of the coordinate system as aligned with the axis of 
the CNT, and the x coordinate is measured along the 
circumference of the nanotube as indicated in Fig. [2^,. 

The four tight-binding wave functions, corresponding 
to a fourfold degenerate single-particle energy level of the 
QD have the forirPSl 



(1) 
where s e (t,l) = (+,-) and v e (K,K r ) = (+,-) 
are spin and valley quantum numbers. Furthermore, a S 
{A, B} is the sublattice index, I is the unit-cell index, 
il cc ii is the unit-cell area, vi a = (xi a ,zi a ) is the position 
of the carbon atom on sublattice a in the Ith unit cell, the 
phase factors-^ in the exponential are <Pk,A — <Pk',b = 
0, <PK',A = Vi an -d fK.B = V ~ n /3 with r\ being the 
chiral angle of the CNT\ and X + = (1, 0) and X- = (0, 1) 
are the two possible spin states with axial polarization. 
The four smoothly varying envelope functions "F^ can 
be obtained by solving the Dirac-type envelope function 



equations^ for v G (+, — ), 



[vf(v x Px + V<TyP z ) + Vconf{z)] 
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Here a x and a y are Pauli matrices, corresponding to the 
sublattice degree of freedom and V con { (z) is a smooth con- 
finement potential, e.g., induced by electrostatic gates. 
Note that our choice of the coordinate system [see Fig. 
2[a)] implies that p z (and not p y ) appears in the envelope 
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Hamiltonian. The functions "F^ and ip vs are normalized, 
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where R is the radius of the nanotube. 

Our goal is to set up a 4 x 4 effective Hamiltonian 
describing the valley mixing due to short-range disor- 
der. Short-range disorder can be caused by any kind 
of atomic faults of the crystalline structure: substi- 
tutional or interstitial atoms, vacancies, adatoms, etc. 
We take into account short-range disorder in the tight- 
binding model as a static random on-site potential Vi, 



(H, 



dis.TB 



ViSij. 



[i = (la) is an index combining 



the unit cell index I and the sublattice index a.] With- 
out the loss of generality, we can assume that the disorder 
potential has zero mean, (Vi) — 0. The short-range im- 
purities are typically charge neutral, and therefore the 
interaction between them is weak. This suggests that 
the random on-site potential is spatially uncorrelated, 
(ViVj) = 8ij(V?). A further plausible assumption is that 
the CNT is homogeneous. Motivated by these observa- 
tions, we model the disorder potential on the different 
sites as independent and identically distributed random 
variables. Since we focus on valley effects, we neglect 
possible sources of spin-dependent short-range disorder, 
such as hyperfine interaction due to 13 C atoms^l and 
adatom-enhanced spin-orbit interactionpB for example. 

To derive an effective 4x4 Hamiltonian describing 
the effect of the short-range disorder, we project the 
tight-binding disorder Hamiltonian -ffdis,TB onto the four- 
dimensional subspace of interest. The corresponding pro- 
jector is 



P = 'Y]\tpva)(lpva 



(4) 



The same method has been used recently by us to analyze 
the effect of hyperfine interaction in carbon-based QDs™ 
The obtained effective Hamiltonian is 



H, 



dis. off 



PH dis , TB P 

(b f + b x T x + b y f y + b z f z ) (g) so 

(6(vfb + & • t) ® so, (5) 



where 



bk — ficell 2^ Vi a F l(J 

la 



(0) 



(6) 



for k G {0.x, y,«}. Here, F^ = EJ*^(n<r)| 2 /2, 
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E^I*r(n.)| 2 /2, and F; o 

The operators fo,x,!/,z are natural representations of the 
Pauli matrices on the two-dimensional Hilbert space 
spanned by ipK and ipK', i.e.) 

fo = IVjtX^jcI + IV'AT'X^JC'I, (7a) 

f x = \iI>k)(iI>k>\ + \iI>k>)(iI>k\, (7b) 

T y = -i\lpK){^K'\+i\^K'){^K\, (7c) 

T Z = \i> K )(ll>K\-\lpK>)(lpK>\, (7d) 

and so is the unit matrix in spin space. The first term in 
Eq. pj, proportional to fo, is just a simultaneous shift 
of the four energy levels. The second term in Eq. (|5| is 
reminiscent of a Zccman coupling Hamiltonian but here 
the roles of spin and the magnetic field are played by 
the valley operator f and the disorder-induced effective 
magnetic field b, respectively. 

In the presence of time-reversal (T) symmetry b z = 0, 
because in this case Tip v is tp- v up to a phase fac- 
tor, implying \^y (r; CT )| 2 = |\ti (rjo-)l 2 and therefore 

(z) 

F} ■ = 0. In CNTs, in the case of moderate magnetic 
fields this statement still holds: an axial magnetic field 
induces an Aharonov-Bohm phase which does not modify 
the electronic wave functions (although induces energy 
shifts), and a perpendicular magnetic field interacts pri- 
marily with the spin of electrons via the Zeeman effect 
as long as the magnetic length is much larger than the 
nanotube radius. Therefore, throughout this paper we 
omit the b z term from the effective disorder Hamiltonian 

Hdis,cS- 

The remaining valley-Zeeman field b = (b x ,b v ,0) is 
random in the sense that it depends on the actual disor- 
der configuration. We consider some statistical proper- 
ties of the valley-Zeeman field in the following. First, the 
averages of its components are 



(b x , y )=n cel lJ2i V l<r) F l 



0,2/) 
la 



= 0. 



(8) 



la 



Second, the variance of the components can be evalu- 
ated assuming that the envelope functions are "flat", i.e., 
1 / (tt cc \\N) , with N being the total number of 
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carbon atoms in the QD. The result is (b 2 . 
thermore, it can be proven that for large N, the quanti- 
ties bo, b x , and b y become statistically independent and 
characterized by Gaussian distributions, potentially fa- 
cilitating future modeling of carbon-based QDs where 
disorder-averaging is necessary. 

Our results imply that the disorder-induced valley 
splittings in the quantum-dot energy spectrum should 



have an order of magnitude \ {b% y ). The typical on- 
site energy Vi on an impurity site is presumably on the 
atomic energy scale, therefore we take 1 eV as an esti- 
mate. Taking a quantum dot containing 10° atomic sites 
and 50 impurity sites we find the disorder-induced valley 



mixing energy scale \/{b 2 y ) « 50/^eV, consistent with re- 
cent experime nts ca rried out on CNT single and double 
quantum dotsP^LLJ 

It is important to note that the model presented in 
this section relies on the assumption that the disorder- 
induced valley-mixing energy scale is much smaller than 
the level spacing in the QD, i.e., the energy distance be- 
tween the fourfold-degenerate level under consideration 
and its neighboring fourfold-degenerate levels. Recent 
measurement d 19 ! 21 ! imply that this assumption is reason- 
able in clean CNT QD devices. In the case when this as- 
sumption is invalid, i.e., if the valley-mixing energy scale 
becomes comparable to the level spacing, then valley mix- 
ing might become efficient between subsequent fourfold- 
degenerate levels, which implies that even the picture 
of independent fourfold-degenerate levels breaks down, 
let alone our model based on the concept of the valley- 
Zeeman field acting on a single fourfold-degenerate level. 
The mechanism that strong disorder mixes subsequent 
levels might actually be a reason for observing twofold 
(as opposed to fourfold) electron shell filling patterns in 
a number of CNT QD experiments. 

For sake of completeness, we give the disorder- 
independent part of the single-electron Hamiltonian cor- 
responding to a fourfold-degenerate QD energy level. 
The finite curvature of the CNT enhances spin-orbit in- 
teraction and leads to a signifi cant splitting (Aso ~ 
100/ieV) of the four levels P^"^ Furthermore, an exter- 
nal magnetic field induces a Zeeman splitting of the spin 
states, and its axial component induces a splitting of the 
valley states as well. Therefore, the disorder-independent 
part of the Hamiltonian is 



„ A so . 

#0 = -7^ T zSz 



i< bB ■ [ ^g s s + ^g v r z z ) . («)) 



where Aso describes the energy splitting caused by the 
curvature-enhanced spin-orbit interaction, g s and g v are 
the spin and valley g factors, respectively, and z is the 
unit vector in the z direction. The form of the Hamilto- 
nian Ho reflects the fact that at zero magnetic field, the 
Kramers theorem implies that the state pairs connected 
by time reversal, i.e., (ipKt,ipK'i) and {ipKi,ipK't) are 
degenerate. In Fig. 2[b), we give an example for the 
evolution of a fourfold-degenerate level with magnetic 
field, which we obtain by diagonalizing H + -ffdis.cff [cf. 
Eqs. ^ and ([9j)]. The parameters used for Fig. [5Jb) 

areP A so = 370/ieV, Jb 2 x + b 2 y = 30/ieV, g s = 2, and 

g v = 54. The valley-independent term proportional to 
&0) which would shift the four levels simultaneously, is 
neglected. 

To conclude this section: we have demonstrated that 
short-range disorder in CNT QDs appears as a random 



(in magnitude and direction) valley-Zeeman field in the 
effective Hamiltonian describing a fourfold- (spin and val- 
ley) degenerate quantum-dot level. We note that our 
derivation is not by any means specific to the particular 
geometry of nanotubes, and we expect the same qualita- 
tive consequences of short-range disorder in the case of 
electrostatically defined QDs in graphene^or siliconPSHSSl 



III. TRANSPORT MODEL FOR A DOUBLE 
QUANTUM DOT 

Our aim in this section is to provide a model for elec- 
tronic transport through a few-electron CNT DQD which 
takes into account the following characteristic features of 
CNTs: (i) fourfold (spin and valley) degeneracy of the 
spectrum (ii) spin-orbit interaction, and (iii) disorder- 
induced valley mixing. In the subsequent sections, we 
use this model to calculate the leakage current through 
a CNT DQD in the spin-valley blockade regime. 



A. Hamiltonian 

We use a constant-interaction Hamiltonian to model 
the few-electron CNT DQD. We take into account a sin- 
gle fourfold (spin and valley) energy level in each QD. 
We consider the case of spin- and valley-conserving in- 
tcrdot tunneling. We write the Hamiltonian in terms of 
creation d\ va (d Rvs ) and annihilation d^ vs (d Rvs ) opera- 
tors of electrons on the left (right) dot having valley and 
spin quantum numbers v and s, respectively, 



y , 



dlA 



dv' 



S,l. 



2-^v ,. 



5i,ss' ^ dvs ^dvs 



and 



ffi 



DQD 



— Hpot + H e 



H so + Hdis + H ma „ n + H u 



H, 



H, 



pot 2^ Crf^d, 

d=L,R 

U 

2 



-ffso — — 



TjT E n d( n d - !) + u 'n L n R , 

d=L,R 

Aso 



E 



Sd,z7~d,z i 



d=L,R 

Hdis = 2^, (bd,xT~d,x + bd,yTd,y) , 
d=L,R. 

-ffmagn = HbB ■ ^ ( -g s S d + -g v 
d=L,R ^ 

#tun = t^2d) Lvs d Rvs +M.c. 



Td.zZ 



(10a) 
(10b) 

(10c) 

(lOd) 
(lOe) 

(lOf) 
(10g) 



The terms (10b)-(10g) in the Hamiltonian describe the 



effects of electrostatic potential difference between the 
dots, electron-electron interaction, spin-orbit coupling, 
short-range disorder, external magnetic field, and in- 
terdot tunneling, respectively. Here d — L, R is the 
QD index, n dvs = d\ vs ddvs, n d = J2vs n dvs, Td,i = 



both Ti and Si (i = x,y,z) are the three Pauli matrices. 
Note that we have incorporated the valley-independent 
disorder-induced terms ~ 17^0 and ~ tr^ into ez, and €r, 
respectively. The operators r d ,k (d = L,R, k = 0,x,y,z) 
defined above are the many-body generalizations of the 
single-particle operator ffe defined in Eq. ph but in Eq. 
( 10 1 and henceforth we suppress the tilde for simplicity. 



We emphasize that the disorder-induced valley- Zeeman- 
fields b L = (b L ,x,bL,y,0) and b R = (b R , x ,b RyV ,0) are dif- 
ferent in general, since the electrons on the left and right 
dot interact with a different set of impurities and there- 
fore feel different disorder configurations. This feature is 
reminiscent of hyperfine interaction in conventional semi- 
conductor DQDs, and will play a critical role in all the 
results we present in the forthcoming sections. 

The system is in the spin-valley blockade regime if 
the available charge configurations for transport are the 
(1,1), (0,2), and (0,1) configurations (in general, these 
numbers might refer to the occupations in addition to 
completely filled shells, see below). In this situation, the 
only relevant parameter of iJ po t + H c ^ c is the energy 
difference (or "detuning") A of (1,1) and (0,2) states: 
A = €l — e R + U' — U. All results presented in this 
work correspond to zero detuning, A = 0, implying that 
the actual values of e^, e R , U, and U' do not affect our 
results. However, since in the following we neglect hy- 
bridization with (1,0) and (2,0), we implicitly assume 
that t/(U - U') = t/(e L - e R ) < 1. 

Motivated by the experiments we try to model here, 
we consider the case of an axial magnetic field: B = 
(0, 0, B). In the CNT QD studied by Kuemmeth et dW 
a spin g factor g s f=a 2 and a valley g factor g v ~ 50 
have been found. However, the results we present in this 
work are insensitive to these values, because (at least 
in the parameter regimes under consideration here) (i) 
spin-Zeeman splitting do not affect the dynamics and (ii) 
we plot the magnetotransport curves against the field- 
induced valley splitting A^ = g v fi R B and not against 
the field B itself. 

Our constant-interaction approximation has the ad- 
vantage of simplicity but also has the drawback that 
it does not account for the recently predicted Wigner- 
molecule formation effect.™^ This restricts the applica- 
bility of our model to (i) short quantum dots, where the 
confinement energy exceeds the interaction energy, or (ii) 
DQD systems where the environment (the metallic gate 
electrodes or the dielectric substrate, for example) pro- 
vides a strong electrostatic screening and hence weakens 
the electron-electron interaction. Wigner-molecule for- 
mation implies a strong suppression of the supersinglet- 
supertriplet gap in the (0,2) charge configuration, which 
suppresses the Pauli blockade as well. The fact that 
Pauli blockade has been observecPHMH in CNT DQDs 
indicates that the samples used in those experiments 
are closer to the constant-interaction regime than to the 
Wigner-molecule regime, which is a further motivation 
for us to use the constant-interaction model in our cal- 



culations. 

Our Hamiltonian does not contain hyperfine interac- 
tion and spin- or valley-flip interdot t unneli ng matrix 
elements, although hyperfine i nterac tio n 1 43 ! 44 ! and spin- 
orbit-induced spin-flip tunneling 45 have proven to be 
important in the understanding of spin blockade ex- 
periments in conventional semiconductor quantum dots. 
We neglect hyperfine interaction in this work because 
theoretical estimates indicate that its characteristic en- 
ergy sca le is b elow 5neV even for fully 13 C-enriched 
samples) j 6 l 10 l 47 l being small compared to other relevant 
energy scales in our system i.e., spin-orbit splitting > 
100/ieV, disorder > lO^eV, interdot tunneling, and val- 
ley splitting > pe V, see forthcoming sections. Note that 
recent experiment d 21 * 22 -! indicate a two orders of magni- 
tude larger hyperfine energy scale than the theoretical es- 
timates, and therefore we cannot be conclusive about the 
relevance of this effect. Experiments in GaAs and CNT 
DQDs have shown that hyperfine coupling becomes espe- 
cially relevant at suppressed interdot tunneling or large 
detuning, therefore our model excluding this mechanism 
might not be adequate in that regime. 

Spin-orbit-induced spin-flip interdot tunneling could in 
principle be present in our system, but only between 
states having a ±1 difference in their circumferential 
quantum number, as it can be deduced from Eqs. (31)- 
(33) of Ref© This possibility is not ruled out in some of 
the spin- valley blockade transport experiments^! as those 
were not performed in the actual (l,l)-(0,2)-(0,l) regime 
(which would imply that in both dots the electrons oc- 
cupy the lowest-energy circumferential mode of the CNT 
and therefore spin-flip tunneling is forbidden). For exam- 
ple, in the (4n+l, 4m+l)-(4n, 4m+2)-(4n, 4to+1) regime 
the spin-valley blockade could take place "on top of" n 
(to) filled shells in the left (right) dot. However, we post- 
pone the analysis of spin-flip tunneling for future work 
and in Sec. |IV| we demonstrate that agreement with ex- 
perimental results can be obtained from our model even 
though spin- flip tunneling is not taken into account. 



to and from the DQD, characterized by the rates Tl 
and Tfl, respectively, ft has the following form: 



Dp = T L J2 [d\ vs pd L vs ~ ^pd L vsd\ vs - ^d Lvs d\ vsP 

VS ^ 

+ T rY1 i d Rvspd\ ivs - -pd Rvs d Rvs - -d Rvs d Rvs p\ 



(12) 



Here the creation and annihilation operators are re- 
stricted to the charge configurations participating in the 
transport process. 



C. Secular approximation 

We assume that the splittings between the eigenval- 
ues of Hdqd are larger than the level broadenings set by 
the tunneling energies hTi, and hT R . This allows us to 
use the so-called secular approximation} 5 ^ i.e., to assume 
that the steady-state density matrix is diagonal in the 
eigenbasis of -ff D Q D . Hence the generalized master equa- 
tion simplifies to a steady-state classical master equation 
(CME), 

= p a = -p a T R y^Pja + T L y^,pjr aj , (13a) 

3 3 

= pi = -pjT L y^ rpj + T R y^ pppip. (13b) 
P P 

Here a, (5 = 1, . . . , 22 (i,j = 1,2,3, 4) refers to the two- 
electron (single-electron) DQD energy eigenstates, p a — 
p aa and pi = pa, and 



r a % = ^2 Ki-ld-Lvsla 

VS 

Via = J^ \(i\dR VS \a 



(14a) 
(14b) 



B. Generalized master equation 

We apply the master-equation formalism to describe 
the transport process through the serially coupled DQD 
system. The DQD charge configurations which are rel- 
evant for the transport process considered here are the 
(1,1), (0,2), and (0,1) configurations. Hence the state of 
the DQD system is described by the 26 x 26 density ma- 
trix p, where the Hilbert space is spanned by 16 states 
in the (1,1) charge configuration, six states in the (0,2) 
charge configuration, and four states in the (0,1) config- 
uration. The time dependence of p is governed by the 
generalized master equation or Lindblad equation 



P= -^[H D qb,p\ + Dp. 



(11) 



The dissipative term Dp describes the tunneling events 



D. Eliminating (0,1) states from the classical 
master equation 

The Hamiltonian has a block-diagonal structure: the 
two-electron [(1,1) and (0,2)] and single-electron (0,1) 
blocks are uncoupled. However, the Lindblad terms do 
couple these sectors because they describe single-electron 
tunneling onto and from the DQD. The coupling is ap- 
pearing in the CME in the form of the rates T^r ai and 

?RPia- 

Throughout this analysis, we consider the case 1\ ^$> 
T R . The reason is that in the spin- valley blockade regime 
the characteristic scale of the rates Ti,r a i are largely in- 
dependent of the spin and valley physics inside the DQD, 
whereas the rates T R p ia are sensitive to those, so in order 
to have the transport via the DQD sensitive to spin and 



valley effects, the outgoing rates TnPi a should provide 
the transport bottleneck. 

We claim that in this limit Tr/T^ — > 0, the steady- 
state CME is reduced to a homogeneous linear set of 
equations Mp = for the vector p = (p%,p2, ■ ■ ■ , P22) 
which contains the diagonal elements of the two-electron 
sector of the DQD density matrix p, and the normaliza- 
tion condition ^2 a p a — 1- The coefficient matrix M is 
given as 






bafiPja 



(15) 



The proof of this statement is a straightforward calcula- 
tion starting from the steady-state CME in Eq. (13 1. 



Having the steady-state occupation probabilities p a 
and the corresponding energy eigenstates |a) at hand, 
we calculate the current as the average decay rate of the 
two-electron states with respect to the steady-state dis- 
tribution, 



1 = eT R y^p a y^p ja . 



(16) 
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FIG. 3: (Color online) Numerical results for the current as a 
function of magnetic-field-induced valley splitting for differ- 
ent values of interdot tunneling (shown). Further parameters: 
A so = 250^eV, b Lx = 20/ieV, b Ly = 10/ieV, b Rx = 80/xeV, 
bn y = OfieV, and therefore b 2 _ /Aso = 23.6 . P oints: numer- 
ical data. Curves: the analytical formula ( 23 1 fitted to the 



numerical data with n* as the single fitting parameter. Lower 
two data sets are scaled as shown. 



IV. STRONG SPIN-ORBIT COUPLING 

In this section, we describe the spin-valley blockade 
effect in a CNT DQD in the case when spin-orbit cou- 
pling dominates the energy spectrum over disorder, inter- 
dot tunneling, and magnetic-field-induced spin and val- 
ley splitting, i.e., Aso ^> b,t,A v ,A s . Here b denotes 
the typical energy scale of the disorder-induced valley- 
Zccman fields on the two dots. The main result of this 

b 2 _ 

section is that we identify a parameter regime (t < ^ — , 
where b 2 _ =b 2 L — b 2 R ) where the current as the function of 
magnetic field (the "magnetotransport curve") shows a 
dip around zero field, and the width of the dip is control- 
lable by the interdot tunneling amplitude t. This field- 
induced increase in th e curr ent is in qualitative agree- 
ment with experiments! 21 * 22 ^ We interpret this result us- 
ing Lowdin perturbation theoryp^ and provide an ana- 
lytical formula for the current which can be well fitted 
to the numerical results using a single fitting parameter, 
the average number of transmitted electrons between two 
blocking events.^ In the following, we describe the case 

b 2 _ , m 

we argue that the findings of 

as 



t 



In Appendix 



A 



this regime can be extended to the regime i« t 
well, and in Appendix |B| we show that they do not hold 

We start our analysis by presenting the numerical re- 
sults for this regime. In Fig. [31 we show the current 
as a function of the magnetic-field-induced valley split- 
ting A„ , for a fixed value of spin-orbit coupling Aso and 
disorder-induced valley fields b^ and b^ (see caption) but 
different values of interdot tunneli ng t. All parameters 
have a realistic order of magnitude! 19 * 21 ! (Note that the 



Zeeman spin splitting A s plays no role in the transport 
process, see below .) In qualitative agreement with recent 
experiments! 21 * 22 * the data in Fig. |3] shows a zero-field 
dip in the current, and the width of the dip is controlled 
by the interdot tunneling t. In all the three cases dis- 
played, the ratio of the zero-field current I = I(B = 0) 



and the maximal current I n 



is 



■fmax/4 ~ 1-5- This 



ratio agrees well with that observed experimentally in 
Ref . I2T1 [see Fig. 3(a) therein], however, in Ref.[22]a ratio 
of Imax/Io ~ 50 has been found [see Fig. 3(e) therein]. 
Below we argue that the factor I msx /Io ~ 1.5 we deduce 
from Fig. [3] is a rough upper bound for this quantity 
in the parameter regime under consideration, and there- 
fore we conclude that our results (i) agree very well with 
the measurement of Ref. [3TJ and (ii) match the mea- 
surement of Ref. [22] only qualitatively, which might be 
due to mechanisms missing from our model or sample 
parameters in the experiment not fitting into the param- 
eter regime we consider here. Further discussion on this 
discrepancy with Ref. [22] is provided in Sec. |VI| In the 
remaining part of this section, we provide an interpreta- 
tion of the numerical results shown in Fig. [3] and derive 
an analytical formula for the current using Lowdin per- 
turbation theory. 

The transition rates in the classical master equation 
[Eq. ( 13 )] are determined by the eigenstates of the two- 
electron Hamiltonian. To provide an interpretation of 
the numerical results, we will describe those energy eigen- 
states using perturbation theory. We start with the two- 
electron Fock basis based on the single-particle states 
ipKfi TpK'i and ipK\,, ipK't [the pairs are energetically sep- 
arated by the spin-orbit energy Aso at zero field, see Eq. 
(|9J|]. The (1,1) states are denoted in the form \Kt,K'f) 
whereas the (0,2) states in the form |0, K t K' f). We 
perform a basis transformation in order to obtain basis 



Spin-orbit energy 
Aso 



Up-spin (S z = +1) 
\K'f,K'f) 



Down-spin (S z = — 1) 
\Kl,Kl) 



Mixed spin (S z = 0) 

±(\Kl,K't)±\K't,Kl)) 

\0,KjK't) 



-^(\Kt,K't)±\K't,Kt)) 
|0,A't K'f) 



^(\Kl,K'l)±\K'l,Kl)) 

\0,KlK'l) 



\0,KfKl) 
-^QK't,K'i)±\K'i,K't)) 

\0,K'fK'i) 



-A £ 



\Kt,Kf} 



\K'i,K'i) 



A s (\Kt,K'l)±\K'l,Kt)) 



■J~2 



\0,KtK'i) 



TABLE I: The 22 basis states used for perturbation calculations in the presence of strong spin-orbit coupling. The six 
(1,1) states involving a minus sign are supersinglets, the ten further (1,1) states are supertriplets. The six (0,2) states are 
supersinglets. 



states which are eigenstates of the two-electron spin-orbit 
Hamiltonian [Eq. (10d)] and have well-defined supersin- 
glet or supertriplet character at the same time. This new 
basis is presented in Table [TJ classified according to their 
properties outlined below. This basis will serve as the set 
of unperturbed states in our perturbation calculations. 

An important simplifying observation is that even in 
the presence of spin-orbit coupling and a magnetic field 
parallel to the nanotube axis, the axial component of the 
electron spin S z is conserved. This allows us to separate 
the 22 states of the two-electron basis to three uncoupled 
spin subspaces (see columns in Table [I]): five states which 
are spin polarized with a polarization aligned with the z 
axis (up-spin states), five states which are spin polarized 
with a polarization antialigned with the z axis (down- 
spin states), and 12 states having mixed spin states. As 
the three different spin subspaces shown in the columns of 
Table IT] are not coupled by any terms in the Hamiltonian, 
the Zeeman spin splitting A s plays no role in the trans- 
port process. Besides their spin state, our unperturbed 
states can also be classified according to their spin-orbit 
energy. Five (five) of those have a spin-orbit energy Ago 
(—Ago), and 12 have a vanishing spin-orbit energy (see 
rows in Table ITJ). 

To visualize the matrix elements of the Hamiltonian, 
in Fig. [4] we show the level diagram of the unperturbed 
basis states we introduced in Table UJ The horizontal ar- 
rangement of the states reflects the charge configuration, 
and the vertical arrangement reflects the spin-orbit ener- 
gies. Red/gray lines denote supertriplet states and black 
lines denote supersinglet states. The green/light gray 
(blue/dark gray) arrows correspond to off-diagonal ele- 
ments of the Hamiltonian in this basis, induced by disor- 
der (interdot tunneling). 

In the polarized spin subspaces [Fig. Qa)] the high- 
and low-energy (1,1) states [dashed lines in Fig. |4^a)] 
are coupled to the (0,2) state via disorder and tunneling, 
resulting in a small decay rate at small fields (A„ <C 



A 



so) 



r. ± = r s i t- 2 



^so 



(17) 



where the ± sign refers to the up-spin and down-spin 

t 2 b 2 

subspaces, respectively, T s = j^Tr <C Tr and b a = 

b^ — bft. Here and hereafter the decay rate of a two- 
electron state a is meant to be the sum of the four tran- 
sition rates into the four different (0,1) single-electron 
states, i.e., J2j=iPja- in the up-spin (down-spin) sub- 
space the decay rate increases (decreases) as the mag- 
netic field increases because the magnetic field pushes 
the high- and low-energy states closer to (away from) the 
zero-energy (1,1) and (0,2) states, cf. Eq. (10f|. The de- 



cay rate T s is fourth order in small parameters, therefore 
we call these four states (i.e., the high- and low-energy 
up-spin and down-spin states) "blocked" . 

To describe the energy eigenstates in the spin-polarized 
zero-energy subspace, conventional degenerate perturba- 
tion theory is not applicable since the perturbative hy- 
bridization of the (1,1) supertriplet state with the (1,1) 
supersinglet would include a zero-energy denominator. 
Therefore, we apply Lowdin perturbation theory^ to de- 
rive an effective Hamiltonian for the zero-energy spin- 
polarized subspace. At zero field, we obtain 



Hq, 



( 



\ 



Ta 



\ 

V2t 
^t j 



(18) 



In -£/o,±i the first index refers to the zero-energy subspace 
and ± to the up-spin and down-spin subspaces. The ef- 
fective Hamiltonian Hq.± corresponds to the following 
ordering of the basis states: (1,1) supertriplet, (1,1) su- 
persinglet, (0,2) supersinglet. Remarkably, H$.± is inde- 
pendent of the angle between the two disorder-induced 
valley fields in the double dot. From Eq. (18) and our 

b - 
assumption ~^- 



t, it follows that the three basis states 



are completely mixed, and each of them acquires a decay 



9 



rate ~ Tr. Therefore with respect to the spin-polarized 
subspaces, we conclude that in the regime considered in 
this section, the ten energy eigenstates can be divided 
to a set of four blocked states decaying with slow rates 
T a> ± <C Tr, and six unblocked states which decay ordcrs- 
of-magnitude faster (with rates ~ Tr) than the blocked 
ones. 

Now we extend this analysis to the 12-dimensional 
mixed-spin subspace [Fig. Kb)]. In the high- and low- 
energy mixed spin subspaces, the effective Hamiltonian 
we obtain is (common diagonal elements are omitted) 



/ 



H±.o — 



\ 





2b 2 _ A„ 





2blA„ 



Ag r 



-4AJ 




V2t 



\ 



V2t 



(19) 



b-Aso 

"SO *■"« / 



Here the ordering of states is analogous to that in Eq. 
( 18 ), and the first index of H± refers to the high- or low- 



energy subspace whereas the second index refers to the 
mixed spin subspace. At zero field, the valley splitting 
is A„ = 0, implying that the first basis state, i.e., the 
(1,1) supertriplet is uncoupled from the other two states, 
in particular, from the (0,2) supersinglet. Therefore, at 
zero field the high- or low-energy (1,1) supertriplet state 
[dashed lines in Fig. Kb)] can decay only due to its 
perturbative coupling (via disorder and tunneling) to the 
two zero-energy (0,2) supersinglets. From Lowdin theory, 
we infer that the decay rate of the high- and low-energy 
(1,1) supertriplets due to these processes is 



r s , = -(r s ,+ + r s ,_) ; 



(20) 



and therefore these two states are blocked in the sense 
defined above. However, according to Eq. (19), a finite 



valley splitting A^ induces mixing between the two (1,1) 
states. This mixing provides an additional decay channel 
for the (1,1) supertriplet state with a rate 



r,= 



26iA 2 . 



i 2 (A 2 



so 



4A? 



rr, 



(21) 



inferred using standard perturbation theory in the field- 
induced coupling term. This rate becomes much faster 
than the slow rate r s _o if the field is strong enough to 



ensure A„ > 



fb a 

V2b 2 _ 



In conclusion, we have found that in 



both the high- and low-energy mixed spin subspace the 
total decay rate of the (1,1) supertriplet state r s! o + T c 
changes dramatically as the magnetic field is turned on: 
at zero field these two states are blocked, having the slow 
decay rate r Sj o whereas at finite field their decay rate 
grows with orders of magnitudes. 

In the six-dimensional zero-energy mixed spin sub- 
space, the order of magnitude of the decay rates is not 
influenced by the magnetic field. Each of these states 
decay fast compared to the slow rate T s . This can be 
derived in the same way as shown at the discussion of 



the zero-energy spin-polarized subspace and Eq. (18) 



Using the explicitly calculated decay rates, we can set 
up a semiphenomenological analytical formula for the 
current. To this end, we regard the transport process 
as an alternation of charge transfer "bursts" (subsequent 
occupation of unblocked states) and blocking events (due 
to occupying one of the blocked states). We assume 
that a burst corresponds to a transfer of n* electrons 
on average, i.e., n* is not necessarily integerP§J Since the 
charge bursts happen fast compared to the time spent in 
a blocked state, the average time between two subsequent 
bursts can be estimated as the average of the decay times 
of the six blocked states, i.e., 



' burst 



2r 



s,+ 



2T-i+2(r a)0 +r c )- 



(22) 



As a burst transfers n* electrons on average, the current 
can be expressed as 



I 



Tbv 



(23) 



Equation d23b, together with Eq. (22) and the decay 



rates calculated above, provides an analytical expression 
for the current as a function of the parameters of the 
Hamiltonian and Tr, having a single phenomenological 
parameter n* . We have fitted this analytical result us- 
ing n* as the single fitting parameter to our numerical 
results (Fig. [3]), and we have found n* « 3.2 irrespective 
of the value of tunneling amplitude t. As seen in Fig. [31 
this value of n* gives an excellent agreement between our 
numerical and analytical results in the considered range 
of magnetic field. By repeating the numerical calcula- 
tions and the fittings for various disorder configurations 
we generally find good agreement between numerics and 
analytics. The values we obtain for n* are typically be- 
tween 1.4 and 5.2, indicating that n* is not universal but 
depends on the details of the Hamiltonian. 

Our analytical result for the current enables us to qual- 
itatively explain two characteristic features of the mag- 
netotransport curves shown in Fig. [3j One of those fea- 
tures is the ratio / m ax/^o ~ 1-5. Evaluating the current 
according to Eq. p3| at zero field, we find Iq — en*r s , 
whereas at high field, where T c ^> T s , we can neglect 
(r^o + r.r 1 in Eq. {22} and find J max « en*6T s /4, 
resulting in the ratio I max /Io ~ 1-5, in correspondence 
with our numerical results in Fig. [3] and the experimen- 
tal data of Ref. HH For this estimate we neglected the 
field dependence of the rates T Si ± but taking that into 
account could only lower the ratio / m ax/^o- A second fea- 
ture observed in Fig. [3] is that the width of the zero-field 
dip of the magnetotransport curve depends on the tun- 
neling amplitude t. This is explained by the fact that the 
crossing-over rate in Eq. (21) depends on the tunneling 
amplitude as T c oc 1/i 2 , i.e., the stronger the tunneling, 
the "slower" the crossover of T c as the magnetic field 
increases, and therefore the wider the zero-field dip in 
the magnetotransport data. In Appendix [A] we argue 
that the conclusions drawn in this section for the case 



t 



-» — can be generalized to the regime t <C a — , and 
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FIG. 4: (Color online) Unperturbed two-electron states (lines) and their energies in a CNT DQD at strong spin-orbit coupling 
and zero magnetic field (cf. Table|T|. Horizontal arrangement of the states reflects charge configuration and vertical arrangement 
reflects spin-orbit energies. Red/gray lines: supertriplet states. Dashed red/gray lines: blocked supertriplet states. Black lines: 
supersinglet states. Green/light gray (blue/dark gray) arrows correspond to off-diagonal elements of the Hamiltonian, induced 
by disorder (interdot tunneling), (a) Up-spin and down-spin states. The indicated x2 degeneracy corresponds to the two 
possible spin configurations. Different spin species are uncoupled, (b) Mixed-spin states. 



therefore the range of validity of our results is actually 

Finally we point out a possible generalization of our 
results. In the system under consideration, the spin and 
valley degrees of freedom play a symmetric role in the ab- 
sence of disorder and magnetic field, since the spin-orbit 
Hamiltonian H so ex s z t z is symmetric in spin and val- 
ley, and the interdot tunneling conserves both spin and 
valley. The results of this section show that if disorder 
provides an inhomogeneous valley-Zeeman field (coupled 
to t x and T y ) in the DQD, then the dynamics becomes 
independent of the spin-Zeeman splitting, and the mag- 
netotransport curve shows a dip at zero axial magnetic 
field. These results can be transferred to the case when 
the role of spin and valley are exchanged: in the hy- 
pothetic case of absence of disorder, an inhomogeneous 
spin-Zeeman field, coupled to s x and s y but not to s z , 
e.g., coming from a perpendicular-to-nanotube-axis mag- 
netic field, would imply that the dynamics becomes inde- 
pendent of the valley-Zeeman splitting, and the magneto- 
transport curve would show a dip at zero axial magnetic 
field. 

To conclude this section: solving the transport mas- 
ter equation numerically, we have found that if Ago 3> 

bt 

2 — ^ t and Ago S> A„, then the magnetotransport 
curves show a zero-field dip with a width that is control- 
lable by the interdot tunneling amplitude t, in agreement 
with a recent experiment.^- Using Lowdin perturbation 
theory, we gave an analytical formula for the current and 
based on that, a qualitative interpretation of the fea- 
tures of our numerical results. We emphasize that the 
observed characteristic magnetotransport pattern is due 
to the different disorder-induced effective valley-Zeeman 
fields on the two quantum dots. Our theory predicts 



a typical ratio of the finite-field and zero-field currents 
Imnx/Io ;$ 1.5, which is in line with the experimental re- 
sult of Ref . HU but different from that of Ref. [531 possibly 
due to experimental sample parameters not fitting into 
the parameter range studied here or mechanisms missing 
from our transport model. 



V. STRONG DISORDER 

In this section, we describe the spin-valley blockade 
effect in a CNT DQD in the case of strong disorder, 
weak interdot tunneling and weak spin-orbit coupling 
(6 ^> t, Ago). In recent experiments on clean nanotube 
QDs, the spin-orbit splitting of the fourfold-degenerate 
ground-state energy level has been found significantly 
larger than the valley mixing energy scale. However, in 
nanotubes with stronger impurity contamination (larger 
radius) the disorder (spin-orbit interaction) energy scale 
is expected to increase (decrease), and the regime con- 
sidered in this section might be reached. A further mo- 
tivation to st udy this regime is its possible relevance for 
certain silicorp2H£5| or g ra phene-based quantum dotspJ In 
those material systems, the spin-orbit interaction is ex- 
pected to be smaller than in CNTs but short-range dis- 
order couples valleys for the same reason as it does in 
CNTs. 

As the main result of this section, we show that in 
the parameter regime under consideration, the magneto- 
transport curve shows a zero-field dip or peak depending 
on the disorder configuration. We find that at a given 
value of the external magnetic field, the current is deter- 
mined by three parameters (if interdot tunneling t and 
emptying rate Tr are fixed): the angle #( tot ) between the 

two total valley-Zeeman fields b L ° = bi, + A v z/2 and 
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FIG. 5: (Color online) Numerical (dots) and analytical (lines) 
results for the current as a function of magnetic-field-induced 
valley splitting for two different disorder realizations. Param- 
eters: t = 10/ieV, bLx = cos#x WOfieV, b^ y = sin#x 100/ieV, 
bR X — 200/ieV and bn y = 0/xeV. Circles: 6 = 7r/4. Boxes: 
= 15tt/16. 
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FIG. 6: (Color online) Two-electron states (lines) and their 
energies in a CNT DQD at strong disorder. Arrows denote 
tunnel couplings between (1,1) and (0,2) states. Decay rates 
of (1,1) states are indicated, cf. Eqs. (|25| and |26[) in text. 
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(a) States and tunnel couplings in the five- 
dimensional spin-triplet subspaces. The same plot refers to all 
three spin-triplet subspaces. (b) States and tunnel couplings 
in the seven-dimensional spin-singlet subspace (cf. Table |H| . 



b 



(tot) 



bji + A v z/2 on the two dots, and the lengths of 

these valley- Zeeman fields 6^° and b R ° . Our analysis 
is analogous to the derivation of the spin blockade leakage 
current induced by hypcrfinc interaction in GaAs double 
dotsp^ with the most important difference being that in 
our case spin-independent disorder provides a blockade- 
lifting mechanism via the valley dynamics, whereas in 
conventional spin blockade, the hyperfine interaction af- 
fecting spin dynamics is responsible for lifting the block- 
ade. 

We start our analysis by presenting numerical results in 
Fig. [5J The figure shows the magnetic-field dependence 
of the current for two different disorder realizations. The 
two curves are qualitatively different: one shows a zero- 
field dip whereas the other shows a zero-field peak. In 
the following, using standard perturbation theory and an 
approximative analytical solution of the master equation 
we show that the qualitative difference between the two 
curves is related to the fact that the angle 9 between 
the disorder-induced valley-Zeeman fields on the two dot 
differs for the two disorder realizations. 

We exploit the fact that the full Hamiltonian commutes 
with the axial (z) component of the electron spin, and 
therefore one can identify four uncoupled spin subspaces 
of the 22-dimensional two-electron Hilbert space. It is 
beneficial to choose a classification corresponding to the 
standard two-electron spin-singlet and spin-triplet states 
\S) S , \T + ) S , \T ) S , and |T1) S (the outer lower index refers 
to "spin" ) . Since these four spin subspaces are uncoupled 
from each other, the spin-Zeeman effect has no effect on 
the dynamics and therefore from now on we disregard 
that. 

In the absence of interdot tunneling, the energy eigen- 
states of the single-electron Hamiltonian are trivial: in 
the left dot they are \K L f),\K L ±), \K' L f), and \K' L 1), 
where \Kl) and \K' L ) are defined as the eigenstates of the 



.(tot) 



(tot) 



2x2 matrix b L • r corresponding to the eigenvalue b L 
and —b L ° , respectively. Single-electron energy eigen- 
states of the right dot are constructed accordingly. The 
two-electron energy eigenstates are the standard Fock ba- 
sis states constructed from these single-electron states 
(as long as interdot tunneling is zero). The resulting 22 
energy eigenstates are classified regarding their energy 
eigenvalue and spin state in Table |TT) 

At this point, we make use of the fact that the mag- 
nitudes of the total valley-Zeeman-fields in the two dots 
are typically different (since they have a random contri- 
bution induced by the random arrangement of disorder) , 
and their difference is typically comparable to themselves 



,(tot) 



/.(tot) 
U R 



,(tot) 



^° t} - 



(24) 



Note that this condition can hold only if the magnetic- 
field-induced valley-Zeeman field A„ does not dominate 
over the disorder-induced component, which restricts the 
validity of the following analysis to the range A„ < b. If 
the condition ( p4| holds, then the separations between 
the seven energy levels considered in Table In] are on the 
order of 6, which is much larger than the interdot tun- 
neling t, and therefore we are allowed to treat t as a 
perturbation, and use the states listed in Table [n] as the 
unperturbed states. 

Due to tunneling, the (1,1) states (upper four rows 
in Table In]) perturbatively hybridize with (0,2) states 
(lower three rows in Table In]) and therefore acquire a 
finite decay rate. As shown in Table [TT] and Fig. [6^,, 
in each of the three spin-triplet subspaces there are four 
(1,1) states and they hybridize with a single available 
(0,2) state. Standard perturbation theory and Eq. ( 14b ) 
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Energy 


|S% 


\T+). 


\To)s 


|T->. 


, (tot) , (tot) 


^5(|^Lt,^.l)-|A-L|,A fl t)) 


\K L t,K R t) 


^(\K L t,K R i) + \K L i,K R t)) 


!^lI,^I) 


l(tot) i .(tot) 


^ 2 (\K' L t,K R i)-\K' L i,K R t)) 


\K' L t,K R t) 


±(\K' L t,K R l) + \K' L l,K R t)) 


I^U,^l> 


,(tot) .(tot) 


^ 2 (\K L t,K' R l}-\K L l,K' R V) 


\K L t,K' R t) 


^(\K L t,K' R i) + \K L i,K' R t)) 


1^1,^1) 


-&L 0t) - 6fl 0t) 


^(\K' L t,K R l)-\K f L l,K' R t)) 


\K' L t,K R t) 


^(\K' L t,K' R l) + \K' L l,K' R t)) 


|#U.*«±> 


2&£ ot) 


\0,K R tK a i) 











^(\0,K R tK' R i)-\0,K R iK' R t)) 


\0,K R tK' R f) 


^ (JO, K R t K' R |) + JO, K R | K' R f» 


|0,A fl |A'U> 


-2&£ ot) 


\Q,K R tK' R i) 









TABLE II: The 22 basis states used for perturbation calculations in the presence of strong disorder. Different columns 
correspond to different spin states and each row has a corresponding energy (left column; spin Zeeman energies are neglected). 
The upper four [lower three] rows contain the (1,1) states [(0,2) states]. 



gives two different decay rates, 



Tit 



2t 2 cos 2 (6»( tot )/2) 



( 



(tot) ,(tot) 



L U R 

2t 2 sin 2 (#( tot )/2) 



.(tot) 



- b ry 



(25a) 
(25b) 



The rate Tit O^t) corresponds to the spin-triplet states 
in the second and third (first and fourth) lines of Table [TT] 
In the spin-singlet subspace, there are three (0,2) states 
to hybridize with, although each (1,1) state hybridizes 
only with two (0,2) states (e.g., \K L t,K R |) - \K L 1 
,Kr t) is n °t coupled to |0, K' R \,K' R \) by tunneling). 
Due to hybridization, the four (1,1) spin-singlet states 
acquire two different decay rates, 



T R 



T2S 

Tr 



= 2i l 



cos' 



(0^/2) ; 2sin 2 (0( tot )/2) 



(tot) ,(tot) 
'L u R 



(tot) , (tot) 
'L T u R 



(26a) 



'It 1 



2cos 2 (6»( tot )/2) sin 2 (6»( tot )/2) 



(' 



(tot) .(tot) 
L U R 



(tot) , (tot) 



'R 



(26b) 



The rate T\s (T<2s) corresponds to the spin-singlet states 
in the second and third (first and fourth) lines of Table [TT] 
We emphasize that the valley dynamics of the spin-singlet 
subspace is remarkably different from the valley dynamics 
of the spin-triplet subspaces and the spin dynamics in the 
spin blockade of conventional semiconductor DQDs^l^Al 
in the latter cases there is only a single available (0,2) 
state to hybridize with, whereas in the former case there 
are three of them. 

From now on we are aiming at deriving an analyti- 
cal formula for the current in leading order in the small 
parameter t/b. As the next step toward that we argue 
that the steady-state occupations of the (0,2) states are 



negligible. There are two facts needed to prove this, (i) 
The steady-state current can be separated to contribu- 
tions from single-electron tunneling via the (1,1) states 
and (0,2) states: I = J2 a e(i,i) P* V a + £ae(o,2) P<* r o<- 
For the (1,1) states, decay rates originate from a weak 
hybridization of the (0,2) states, therefore in the first 
sum, T a ~ (t/6) 2 r_R < T R . For the (0,2) states, decay 
rates come from direct coupling to the right lead, hence 
in the second sum, T a ~ T R . (ii) The (0,2) states are 
"difficult to load" and "easy to empty", therefore, as it 
can be shown rigorously, their steady-state occupations 
are ~ (t/b) 4 whereas the occupations of the (1,1) states 
are ~ 1. As a result concerning the current, this means 
that the contributions from the (1,1) states provide the 
leading-order result, and the (0,2) states can be elimi- 
nated from the classical master equation [Eq. (15)] . 

The steady-state CME retrieved after the elimination 
can be solved analytically using the ansatz 

Pa= ^ F " ==1 , [a 6(1,1)], (27) 



y , 



. r: 



'£(1,1) 

which expresses that the occupation probability of a state 
is proportional to the lifetime of that state. The fact that 
this simple ansatz solves our classical master equation is 
a consequence of the equivalence of the 16 (1,1) states in 
the following sense: if one of those is filled randomly with 
a uniform distribution, then after one transport cycle the 
occupations are still uniformly distributed. Mathemati- 
cally, the 16 x 16 "return probability matrix" of the (1,1) 
states 

R aP = ]T " aJ *" , [a,/? 6 (M)] (28) 

36(0,1) ^a' r a'j l^j'Pj'P 

is doubly stochastic. The key observation in proving this 
is that the row sums of the (1,1) — > (0, 1) transition prob- 
ability matrix Pjp/ y^i Pj>p are equal, which is a conse- 
quence of the vanishing detuning between (1,1) and (0,2) 
states. Note that the connection between R and the co- 
efficient matrix M of the CME is 



R a /3 — I 5, 



M, 



a/3 



a/3 



J2ke(o,i)Pk/3 1 



(29) 



a,;8e(l,l) 
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Using the solution in Eq. (27) and the rates listed in 



Eqs. (25 1 and (26 1, we obtain an analytical formula from 



Eq. (16) for the steady-state current through the DQD, 



I = 



16 



6 r 



IT 



IT) 



2 r 



l.S 



r 



(30) 



■is , 



This analytical result is compared to numerical results 
in Fig. [5] and a good correspondence is found. Our an- 
alytical result gives an insight on how the angle 6*( tot ) 
influences the current: at angles close to and 7r, i.e., 
at parallel and antiparallel valley-Zeeman fields on the 
two dots, the current is suppressed, since then either Tit 
or r 2 T is small and that makes the denominator in Eq. 
(30) large. This qualitatively explains the zero-field dip 
in Fig. [5] in the case of 6 = 157r/16: as the magnetic 
field and hence A„ is increased, the angle #( tot ) crosses 
over from 157r/16 toward (since the magnetic field is 



,(tot 



(tot)x 



enforcing alignment of b L and b R ), starting from 
and ending at suppressed current values, but sweeping 
through a region of enhanced current. 

Our conclusion of this section is that in the considered 
regime the magnetotransport curve shows either a zero- 
field dip or peak, depending on the disorder configura- 
tion. Although our analysis in this section was based on 
the complete absence of spin-orbit coupling and detuning 
between (1,1) and (0,2) states, and a perfect alignment of 
the magnetic field and the CNT axis, we expect no quali- 
tative changes in the results in the case of weak spin-orbit 
interaction and detuning Aso , A <§C 6 and/or a small mis- 
alignment of the field, since those factors have no effect 
on the tunneling amplitude and can only slightly mod- 
ify our unperturbed basis states and the corresponding 
energies. 

As discussed in Sec. IV in recent experiment ;! 21 ! 22 ! 
with CNT DQDs a zero-field dip has been found with 
a dip width controllable by the tunneling amplitude. In 
the analytical and numerical results of this section, the 
width of the dip (predicted for certain disorder configu- 
rations) is insensitive to the tunneling amplitude because 
t appears in the current as a t 2 prefactor only. Therefore, 
we conclude that the parameter regime of the measure- 
ment was probably different from the one considered in 
this section. This opinion is supported by the facts that 
in Ref. [5T|the ratio Ago/fr ~ 7 has been estimated and 



that in Sec. IV in a different parameter regime we have 
found a qualitative agreement with experiments. 



VI. CONCLUSIONS 

We have found that in the regime of strong spin-orbit 
interaction, the magnetic-field dependence of the leakage 
current shows a zero-field dip with a width tunable by 
the interdot tunneling amplitude t (Sec. 



IV I. We have 



shown that the ratio of the finite-field and the zero-field 
current is typically I max /Io ^ 1-5- Both the trend in the 
magnetotransport data (i.e., the zero-field dip) and the 



numerical value 1.5 agree well with those found in Ref.l2"T1 
[see Fig. 3(a) therein]. In the measurement of Ref. [22] 
[Fig. 3(e) therein] the qualitative behavior is similar to 
our prediction, but a much larger ratio, I m a,x/Ia ~ 50 has 
been found, which deviates significantly from the predic- 
tion of our model. This deviation might arise from the 
parameters of the measured sample not fitting into the 



parameter regime considered in Section IV (In Ref 



the values of spin-orbit interaction energy and valley- 
mixing energy have not been estimated. Interdot tun- 
neling has been quoted as t ~ 50^eV.) Another potential 
reason for the deviation might be that certain features 
and mechanisms possibly important in the spin-valley 
blockade are excluded from our model. A relevant mech- 
anism might be the spin- and/or valley-relaxation due 
to electr on-phonon i nteraction or electron exchange with 
the leads J51i. i | i3 | 46 | 60 | jf those relaxation rates are compa- 
rable to or larger than the lead-dot tunneling rates then 
they could affect the transport properties. Another pos- 
sibly influential effect disregarded in our model might be 
the emergence of strongly correlated Wigner-molcculc- 
like states due to the strong electron-electron interaction 
in nanotubespH21 which would imply the rearrangement 
of the energy level structure shown in Table [I] and Fig. [4] 
and therefore could lead to a qualitatively different trans- 
port behavior. 

We have studied the influence of disorder on the spin- 
valley blockade in the case of small spin-orbit interaction, 
where the dominant energy scale is that of the short- 
range disorder (Sec. Ivl). In this regime, the leakage cur- 
rent can show a zero-field dip or peak, depending on the 
disorder configuration. Although we are not aware of any 
measurements carried out in this regime, we think that 
our results might be relevant for future experiments on 
graphene- and silicon-based double quantum dots. 

Our QD model incorporating the valley-mixing effect 
due to disorder can serve as a starting point for future 
theoretical work on CNT QDs. The fact that the valley- 
mixing effective Zeeman field depends on the electronic 
wave function and the disorder configuration felt by the 
electron implies that this valley-Zeeman field changes as 
the electron is replaced. This feature might allow for 
resonant electronic valley manipulation similar to recent 
spin manipulation experiments in conventional semicon- 
ductor QDs using spin-orbit coupling^ and hyperfine 
interaction.-^ Note that in clean nanotubes, such reso- 
nant techniques do not require a magnetic field since the 
two valley states having the same spin are split by the 
spin-orbit splitting even at zero field. A further possible 
application of our model could be to describe the pulsed- 
gate experiments of Churchill et a/p2 which intended to 
infer relaxation and decoherence times of two-electron 
spin- valley states in a nanotube double dot. As those 
results have been obtained using isotope-enriched sam- 
ples, incorporating the spin- and valley-mixing hyperfine 
interaction's might also be necessary. 

In conclusion, we have found that spin-independent 
short-range disorder in carbon nanotube double quantum 



dots can lift the spin-valley blockade. In our transport 
model, we account for valley degeneracy, spin-orbit en- 
ergy splitting, and disorder-induced valley mixing, which 
are characteristic features of nanotubc quantum dots dis- 
tinguishing them from their counterparts in conventional 
semiconductors. The main result of this work is that 
in the regime of strong spin-orbit interaction our model 
predicts a zero-field dip in the magnetic-field dependence 
of the leakage current, with a dip width tunable by the 
height of the interdot tunneling barrier. This behavior 
is in accordance with recent experiments. Our analysis 
for the regime of strong disorder, which has possible rel- 
evance for graphene- and silicon-based double quantum 
dots, predicts that the magnetotransport shows either a 
zero-field dip or peak depending on the disorder configu- 
ration. 

Note added in proof. After completion and submission 
of this work we became aware of two related theoretical 
studiespS^l w hich focus on spectral properties of per- 
fectly clean (disorder-free) CNT DQDs. The work of von 
Stecher et alW^ also describes the effect of strong corre- 
lations on transport in the Pauli blockade regime. The 
purpose of the present work is to account for disorder- 
induced effects in the spin-valley blockade, which makes 
it distinct from Refs .loll and l62l 
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there is an important deviation in the procedure com- 



pared to Sec. IV In the zero-energy spin-polarized sub- 
spaces [described by H ± in Eq. (18)], the three basis 



state does not mix evenly but instead the (1,1) super- 
triplet and supersinglet states hybridize strongly with 
each other and these hybridized states themselves hy- 
bridize only weakly with the (0,2) supersinglet via tun- 

neling. This leads to decay rates ~ If Tr <C T^ (in 

contrast to Sec. IIVI where these rates were found to be 
~ r^). An important point is that if b a ~ 6_, which is 
typically true due to the random nature of the disorder- 
induced valley fields b^ and b_R, then the decay rate 

A 2 t 2 

~ — 8? — r# is still orders-of-magnitude larger than the 

decay rate T s of the blocked states, and therefore these 
states can still be considered as unblocked. The same 
argument applies for the six states in the zero-energy 
mixed-spin subspace, therefore each of those can be clas- 
sified as unblocked. 

b 2 _ 

From these results we conclude that in the case -» — 3> 

A so 

t the set of blocked states (which determine the character 
of the magnetotransport curve) and the form of their 

b 2 . 

decay rates are the same as found in the regime ^ — ~ t, 
and hence the conclusions drawn there hold here as well. 
This finding is confirmed by numerical calculations (not 
shown) . 
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Appendix A: Strong spin-orbit coupling and ^ ^> t 

We have shown in Sec. |IV| that in the regime of strong 

b 2 _ 

spin-orbit interaction and ^ — ~ %■> the magnetotrans- 
port curve develops a zero-field dip with a width which 
is controllable by the tunneling amplitude t. Here we 
argue that the above condition can be relaxed and the 

b 2 _ 

same statement is true in the regime ^ — 3> t (provided 

A S o>6,i,A„). 

We recall that in Sec. IIVI we have found that at zero 
field six out of the 22 two-electron states are blocked, and 
two out of those six show a blocked-unblocked crossover 
as the magnetic field is switched on. We have determined 
the field-dependent decay rates of these six states without 

v t, therefore those results 
3> t as well. 



using the assumption 



Aso 
b 2 



hold in the current case 

The effective Hamiltonians in Eqs. (jl8|) and ( |19| are 
also valid in the current case. However, from this point 



Appendix B: Strong spin-orbit coupling and -^- 7 — <C t 
Here we argue that the zero-field magnetotransport 



dip, discussed in Sec. IV and Appendix \A\ gets smeared 
out if the tunneling amplitude is increased to the regime 

.2 

-£=- -C t such that t S> 6_^/A so /6 a - (The latter con- 
dition is stronger than the former one provided that 
t>L 7^ b-R and therefore b a ^ 0.) As before, we restrict 
the discussion to the strong spin-orbit coupling regime: 
A so >M,A. u . 

As pointed out in the analysis after Eq. (19), the 



blocked-unblocked crossover of the high- and low en- 
ergy mixed-spin supertriplets, which gives rise to a zero- 
field dip in the magnetotransport, occurs around the 



magnetic field where A„ = 



t 2 b a 

V2b 2 



This implies that if 



t ^> 6_-\/Ago/&a, then the crossover would take place 
only in the high-field regime A„ 3> Ago and not in the 
low-field regime A„ <C Ago under consideration. As a 
consequence, the character of the low-field magnetotrans- 
port curve will be determined by the field-induced evolu- 
tion of the slow decay rates, resulting in a parabolic peak 
around zero field. This finding is confirmed by numerical 
calculations (not shown). 
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